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ABSTRACT 

Method This invention is a method of separating induced microseismicity signals from 
seismic signals acquired within the context of active seismic monitoring operations 
carried out in underground zones under developmen t which has application for 
underground hydrocarbon or fluid storage reservoir , 

- Seismic records are formed from signals emitted by one or more seismic sources 
controlled by orthogonal signals. In this case, the signals are processed to separate the 
respective contributions of the seismic sources to the signals received and to reconstruct 
the seismograms equivalent to those that would be obtained by actuating the seismic 
sources separately. The induced microseismicity signals (passive seismic) are separated 
from the seismic signals coming from sources & (active seismic) essentially by isolating 
the contribution thereof by comparison with a spectral reference model at the 
fundamental frequencies emitted and at the respective harmonics thereof, and by 
reconstructing in the time domain the induced microseismicity signals and the signals 
coming from seismic sources 5. 

- Applications : seismic monitoring of underground hydrocarbon or fluid storage 
reservoirs. 



BACKGROUND OF THE INVENTION 

FIELD OF THE INVENTIO NF ield of the Invention 

rOOOn The invention relates to a method intended for discrimination of 

discriminating m icroseismic signals induced by the development of an underground 
reservoir, from among seismic signals emitted within the context of seismic monitoring 
operations. 

Description of the Prior Art 

r00021 In order to follow the evolution of an imderground reservoir under 

development, it is well-known to permanently install seismic sources and seismic 
receivers, at the surface or in wells, whether for r eservoir exploration or for 
development wells, and to carry out, at predetermined time intervals, recording of the 
seismic signals from the formation. These signals may be seismic signals induced by the 

site development activity or reflected by the geologic discontinuities of the formation in 

r 

response to the emission of disturbances in the ground by means of one or more seismic 
sources. Comparison of the recordings performed at different times of the site 
development provides information on the evolution thereof 
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BACKGROU>JD OF THE INVENTION 

r00031 ^Various long-term seismic monitoring systems are described for example 

in the following p atents EP-591,037 (US-Eatent_5,461,594), FR -French patent 
2,593,292 (US -patent 4,775.009), PR ^and French patent 2,728.973 (US-5,724,311) or 
FR-2,775,349. 

r00041 P atet^EP -patent 748.457 (US -patent 5.724.3 in filed by two of the 

applicants describes a permanent monitoring system allowing to obtain obtaining 
perfect reproducibility of the operating conditions in repetitive seismic monitoring 
operations carried out in an underground zone traversed by at least one well or borehole, 
notably in an xmderground gas storage reservoir. This system comprises, installed on a 
permanent basis, one or more set(s) of seismic receivers (buried at the surface or placed 
in one or more wells), several repetitive seismic sources (buried or at the surface), and a 
permanent linking network for selective energy supply to these sources. A central 
station remote controls selectively each one of these sources and records the seismic 
signals coming firom the underground zone in response to the seismic waves transmitted 
selectively in the ground by the sources. 

rOOOSI ^AU these permanently installed sources whose coupling with the 

surroxmding formations remains stable and this supply network, at least partly buried 
and whose surface coverage area is limited, allow to carry allows cairving out a series 
of long-term seismic monitoring operations imder stable operating conditions, without 
any risk of incompatibihty with the development site activities. 

rOOOSI P atent FR French patent 2.775.349 (US -patent 6.182.012) filed by the 

same two applicants describes a method and a device intended for permanent 



3 



monitoring of an underground formation in which cavities are formed for seismic 
sources (preferably deep enough to be coupled with the formation below the weathered 
zone (WZ)), at least one well is drilled either below each source or in the immediate 
vicinity thereof An antenna conaisting of h aving a pluraUty of seismic receivers 
(geophones, hydrophone, etc.) is lowered into each well thus drilled and coupled with 
the surrounding formation. With this device, emission-reception cycles are carried out 
with emission of seismic waves in the formation by at least one source and acquisition 
of the signals picked up by the seismic receivers, in response to the waves reflected by 
the formation. The receivers can for example be associated with a fluid transfer tube 
lowered in the well in order to connect an underground zone to a surface development 
device and coupled with the formation surroimding the well. 

r00071 The various sources of the seismic monitoring device can be actuated 

successively with a sufficient time interval between the triggering times for reception of 
the waves reflected by the investigated zone being investigated . It is also possible to use 
sev e ral seismic sources emitting the same signals that are triggered simultaneously to 
increase the power emitted. 

rOOOSI Pat e nt FR French patent 2,805.051 filed by the applicants describes a 

method and a system int e nd e d for seismic monitoring of an underground formation. As 
diagrammatically shown in Figs.l to 4, the system comprises for example a network of 
seismic antennas 2 consi s ting each ef -having a series of seismic pickups 4 arranged at 
regular intervals along a well 3 drilled in the ground. This network can be regular as 
shown in Fig.2 or irregular. The pickups are, for example, vertically oriented one- 
directional or multi-axis (triphones) geophones and/or hydrophones. A seismic soiu^ce 5 
is arranged near to each antenna 2. Piezoelectric type vibrators such as those described 
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in French p atent FR-2,791,780 (US -patent 6338,394) filed in the joint names of the 
applicants are advantageously used as sources and permanently installed in the 
immediate vicinity of each antenna 2 . An electronic control and recording module 6 
(Fig.l) can be associated with each antenna. These electronic modules 6 can be 
connected to a central control and synchronization laboratory 8. According to another 
implementation mode diagrammatically shown in Fig.4, antennas 2 are directly 
connected to a central laboratory 8 including all the individual electronic control and 
recording modules 6 (Fig.l). 

r00091 Seismic waves are emitted in the underground formation by one or more 

seismic sources (vibrators) 5. In the case of multiple sources, the latter emit 
simultaneouslv T and tfeey-are controlled by orthogonal signals so as to form a composite 
vibrational signal. The vibrators work permanently or at least for a sufficiently long 
time to obtain a sufficient signal-to-noise ratio. The seismic waves generated by seismic 
sources 5 are propagated downwards (downgoing waves 9). These incident waves are 
first recorded by receivers 4 of each well 3. The waves reflected by the discontinuities 
of the zone (seismic interfaces) are propagated upwards. These upgoing waves 10 are 
also recorded by the various receivers 4. The upgoing and downgoing waves are thus 
superimposed in the seismograms. The records are subjected to the conventional 
processings known to th e man sidll e d in the art suited for those obtained by means of 
the VSP technique (Vertical Seismic Profiling). Finally, the respective contributions of 
the seismic sources to the composite vibrational signal are discriminated and the 
seismograms equivalent to those that would be obtained by actuating the seismic 
sources separately are reconstructed. Sinusoidal signals of different frequencies, in their 
fundamental components as well as in their respective harmonics, or signals based on 



wavelets, Legendre polynomials or random series, etc., are for example used as 
orthogonal signals. 

[00101 In the case notably where the orthogonal signals emitted are sinusoids, 

discrimination of the respective contributions of the seismic sources is for example 
performed by determining the amplitude and the phase of the composite vibrational 
signal at the fundamental frequencies of the control signals applied to the seismic 
sources. This discrimination comprises for example weighting of the signals recorded 
by a bell weighting (or apodizing) factor and determination of the amplitude and of the 
phase of the composite signal, or selection by Fourier transform of lines of the complex 
spectrum associated respectively with the various weighted signals. Reconstmction of 
the seismograms corresponding specifically to the various seismic sources is carried out 
for example by applying, after separation thereof, an inverse Fourier transform to the 
lines associated respectively with the various weighted signals. 

room ^Furthermore, it is well-known that reservoir development induces 

mechanical stress variations in the reservoir and in its overlying bed. These stress 
variations have the effect of generating seismic signals that are propagated in the 
formations. The receivers of the various antennas 2 record these signals (upgoing 
waves 11). 

r00121 Since seismic sources 5 are arranged much closer to receivers 4 of 

antennas 2 than to the reservoir and since the energy they emit that is emitted is much 
more intense than the energy of the induced microseismic signals generated in the 
reservoir, the microseismic signals are generally concealed and undetectable. 



SUMMARY OF THE INVENTION 

r00131 ^The method according to the invention allows to s e parate separation of 

induced microseismicity signals and multisource seismic signals within the context of 
monitoring operations in an underground zone under development. 

r00141 It comprises emission of seismic waves in the formation by coupling 

therewith one or more seismic sources. In the latter case, these sources emit 
simultaneously orthogonal signals so as to form a composite vibrational signal, 
reception of the signals reflected by the formation in response to the emission of seismic 
waves, recording of the signals received by at least one seismic pickup and processing 
of the recorded signals record e d so as to separate the respective contributions of the 
seismic sources to the signals received and to reconstruct the seismograms equivalent to 
those that would be obtained by actuating the seismic sources separately. 

r00151 The induced microseismicity signals (passive seismic) are separated from 

the seismic signals emitted essentially by isolating the contribution of each source by 
comparison with a reference spectral model at the emitted fundamental frequencies 
emitted and at the respective harmonics thereof, and by reconstructing in the time 
domain the induced microseismicity signals and the seismic signals emitted. 

r00161 The spectral contribution of the microseismicity signals to the spectrum 

of the signals received is obtained for example by subtracting the amplitude and phase 
values associated with the reference spectral model from the amplitude and the phase of 
the spectrum associated with the records. 
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r00171 According to an implementation mode, a cxurent model fomied by 

updating a spectral model taking account of the contribution of previous recording 
cycles is for example used as the reference spectral model. 

r00181 ^The current spectral model can be formed by determining a mean value 

of the frequency spectra formed from (earlier and/or later) records obtained for the same 
source and the same frequencies. 

r00191 ^According to another implementation mode, the current spectral model is 

formed by determining a median value of the frequency spectra formed from earlier 
records obtained for the same source and the same frequencies. 

r00201 According to another implementation mode, the current spectral model is 

formed by extrapolation or interpolation from the frequency spectrum, from close 
spectral values. 

r002n ^According to another implementation mode, the current spectral model is 

formed by using the line obtained at this frequency. 

r00221 ^According to an implementation mode, discrimination of induced 

microseismicity signals from among signals emitted within the context of active seismic 
monitoring operations is obtained by means of the following stages : 

a) for each recording n of a recording cycle p, the respective contributions of the 
various sources at the fiindamental frequencies are calculated^! 

b) the ratio of the contribution to a current spectral model formed by updating a 
previous spectral model from frequencies emitted during the previous recording and 
from their harmonics is then calculatedr: 
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c) the part of the recording n of cycle p that can be associated with the active 
seismic monitoring operations is deduced^; 

d) the part of the recording n of cycle p strictly relevant to the passive microseismic 
activity is deduced?; 

e) the seismograms that can be associated with the active seismic monitoring 
operations are formed by inversion in the time domain of the respective spectral 
contributions of each seismic source at the fundamental frequencies and at their 
harmonics?; and 

f) the underlying microseismic signals contained in the records are formed by 
inversion in the time domain from the part relevant to the passive microseismic activity. 

r00231 The spectral contribution is obtained by multiplying a transfer fiinction 

between a wavelet characteristic of the source and a seismogram associated with a given 
receiver by a wavelet characteristic of the source. This transfer ftmction can be 
continuously updated during a current cycle from an estimation made during a previous 
cycle and from an initial estimation made during the current cycle by the relation. 

BRIEF DESCRlPTIO>[ OF THE FIGURES BRIEF DESCRIPTION OF THE 

DRAWINGS 

r00241 Other features and advantages of the method according to the invention 

will be clear from reading the description hereafter of an embodiment given by way of 
non limitative example, with reference to the accompanying figures wherein : 

r00251 Figxu-e F ig. 1 diagrammatically shows the monitoring device installed in 

the field for seismic monitoring of an underground formation^; 
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r00261 Fipur e Fig. 2 shows a mode of layout of the emission-reception units in 

the fieldj; 

f00271 Figur e Fig. 3 diagrammatically shows one of these emission-reception 

units comprising a seismic source and a reception antennar; 

f00281 Figui e Fig. 4 shows a variant of the monitoring device of Fig. and 

[00291 — Figur e Fig. 5 shows examples of signals and frequency spectra 

mentioned in the description. 

DETAILED DESCRIP TION DETAILED DESCRJPTION OF THE INVENTION 

[00301 ^The method tbus- of the invention allows to carry carrying out seismic 

monitoring operations in an underground zone by means of a series of seismic pickups 
and one or more vibrational seismic sources. If several simultaneous sources are used, 
they are actuated simultaneously by signals at different frequencies selected so as to 
allow discrimination of the contributions of each source in the seismograms obtained 
from the signals received and recorded. 

[00311 This is generally performed by control of the various sources by 

sinusoids of different frequencies and by means of known numerical calculation 
techniques such as inverse Fourier transform to separate the contributions to the 
seismograms obtained from the various seismic sources. 

[00321 ^Distortion caimot be disregarded with mechanical sources. 

Simuhaneously with a frequency fi, a single source Sj emits frequencies 2fi, 3fi ... nfj. 
Consequently, if fi and fj are the respective frequencies of the two sources Si and Sj of 
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the series of sources, we must have ^ % as well as fi 2fj, fj ^ Sfj, ... ft nfj in order 
to be able to separate their respective contributions. 

r00331 As already described in the aforementioned French p atent FR-2,805,051, 

if the composite signal Pe consisting of the sum of N sinusoids {fj, Ai, with 1 < i < N 
is emitted, all the frequencies being contained in a spectral band contained between two 
boundary frequencies fb and ff, the seismogram Tt observed at point R will have as 
Foxuier transform at the frequency fu the number of module Aj and of phase Oi equal to 
the amplitude and to the phase of sinusoid Tj. It is thus possible to reconstruct by 
inverse Fourier transform seismogram Tt by emitting successively all the sinusoids of 
frequencies fb to ff. 

r00^41 In the case where, for example, all the amplitudes Ai are 1 and all the 

phases Oj are zero, the signal Pt obtained is very close to the signal resulting from 
autocorrelation of a sliding-frequency signal contained in the sweep range [fb — fJ, 
commonly used in vibroseismic prospecting. According to the discrete Fourier 
transform theory well-known to the man skill e d in the art, if one wants to listen to 
source Si during the time tc, the frequency increment between the sinusoids is Af = I/tc 
and the number of sinusoids necessary is Nf = (ff - fb)te. 

r00351 N vibrators SI, S2, Si,... Sn installed in the field can thus be 

simultaneously excited by means of vibrational signals whose frequencies are such that 
each source is excited successively by each above sinusoid Nf at any time, on condition 
that the respective frequencies of the sinusoids emitted at the same time by the various 
seismic sources are all different from one another. Nf recordings corresponding to the Nf 
frequencies selected in the sweep range are thus performed with each vibrator. 



Separation of the signals received by the pickups in the field, in response to the 
simultaneous emission of the various signals, is thus obtained by selection of the line at 
the suitable frequency. 

r00361 Cycles of Nf recordings are thus repeated. The record E„ of a given 

acquisition cycle contains a set of fundamental and harmonic frequencies which only 
depends on n. Starting from this observation, the method that will be applied to detect 
the effects of microseismicity is based on the observation which has proved true in 
practice according to which the amplitude and phase variations observed at the 
frequency f and at the time t are very little different from those observed at the 
frequency f - Af and at the time t - At. This allows to form forming a cumulative spectral 
model constructed from records of previous cycles at the same frequencies but fomied 
on a long time diuing which the characteristics of the seismic signals emitted are likely 
to vary, and pemianently adapted from records obtained during a sufficiently short 
period for the variations of the signal to be disregarded. The efficiency of the proposed 
method i s o ptimxrai w hen t he s uccessive f requencies e mitted b y t he s ame s ource a re 
regularly incremented. Cycles in upgoing and downgoing frequencies can for example 
be altemated. 

r00371 ^The p rocedure f or s howing i n p ractice t his s pectral c ontribution o f t he 

underlying microseismic signals is as follows. 

^Notations : 

■ fi.„ designates the fimdamental frequency emitted by source i during the recording 
n (this frequency is independent of the cycle). 



In order to simplify the notations of the following quantities, the subscript giving 
the number of the receiver is left out. Whatever this number, the quantities are 
calculated in the same way. 

■ Rp.„ designates the recording of rank n of cycle p, 

■ Cp.i.„ the contribution of source i to the recording n of cycle p, 

■ Mp,„ the initial model of the active part of the recording n of cycle p. This model 
contains only the frequencies emitted and their harmonics during recording n, 

• Ep.„ the ratio between the active contributions of the recording n of cycle p and the 
initial model of the active part of the recording n of cycle p, 

■ Ap^ the active part of the recording n of cycle p, 

■ Pp,„ the passive part of the recording n of cycle p (microseismic), 

■ k the differences updating coefficient (typically of the order of 15 to 25 %), 

■ h the models updating coefficient (typically of the order of 5 to 10 %). 

r00381 Considering these notations, the method first consists, for each recording 

n of cycle p, in calculating the respective contributions Cp,,.„ of the various sources at 
their fiindamental fi-equencies and at their harmonics by means of the method described 
in the aforementioned Frenclij)atent FR-2,805,051. Owing to the possible presence of 
microseismic energy, the contribution of the various sources is likely to be affected by 
noise. 

[00391 This estimation can be improved by associating the contributions of all 

the receivers of a same recording Cp.,,„. These contributions can be considered to be the 
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product of two factors among which one depends only on the soiirce (and not on the 
receiver) and the other on the receiver, and not on the source. 

r00401 ^If subscript r represents the receiver, 

T only depends on p, i and n by the variation of the characteristics of the receiver 
and of its environment with time. These variations are very slow and can be disregarded 
or estimated more precisely. 

[00411 S and T are determined to within one factor. The average of the 

downgoing waves can thus be selected for Sp.i.„. This average is obtained by means of a 
calculation well-known to tlios e skilled in VSP processing (flattening of the downgoing 
waves and average). The Tpj,„^ are then the transfer functions between the downgoing 
wave and the seismogram. 

r00421 In cases where the variation of T is to be taken into account, updating 

coefficient h can for example be used to update Tp,,,„,, fi-om Tp,i^.i^ as follows : 

r00431 ^If C Op.,^., d esignates t he initial e stimation o f c ontribution C and T Op,,^., 

that of the receiver factor, S is calculated by flattening and averaging of the downgoing 
waves, then TOp.,^,, is calculated by the relation : 



Co . 
To . =— 



[00441 ^The final estimation of the receiver factor is obtained by 

Tp.i,„.r = (1 - h)Tp,,,„.,,r + hTOp,,,„,r 
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and the final contribution of source i to cycle p of recording n by : 

We then caloulat e Then the ratio Ep^ of the current recording to the spectral model 
by the relation is calculated : 

r00451 This calculation requires knowledge of Mp,„ whose calculation described 

below is carried out during cycle p-1. It comprises smoothing oh— the previous 
frequencies so as to limit the effect of the microseismic energy in its estimation. The 
proposed formula only uses the previous records and thus allows real-time conditions. 
In difficult cases, this sum can be replaced by an average or a ciurent median. Deferred- 
time conditions are then necessary (the time required to carry out the second half of the 
recordings used in calculation of the average). 

r00461 ^If Mp+1^ represents the initial model of the active part of the recording n 

of cycle p+1, then : 

r00471 ^This model essentially consists of a weighted weights an average of the 

current contribution and of the previous contributions. The formula used allows 
variation of the model : the contribution of cycle p-q to the model used for cycle p is 
multiplied by coefficient (l-h)** which tends to zero all the faster as h is close to 1. This 
allows model Mp^ to follow the seasonal variations, and ratios Ep^ to reflect only the 
daily variations. Knowing spectral model Mp^ and ratio Ep,a of the contribution to the 
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previously established model, w e d e duce p art Ap,„ of the recording n of cycle p relevant 
to the active seismic part is deduced : 

A = M E 

r00481 Similarly, part Pp.„ of the recording n of cycle p strictly relevant to the 

(passive) microseismic activity is calculated by the relation : 

r00491 ^The respective spectral contributions Cp.i.„ of each source Si at the 

fundamental frequencies and their harmonics, after completion of a cycle of 
measurements, provide, by inverse FFT, the seismograms (in the time domain) used 
within the context of active exploitation. 

rOOSOl The estimation of the underlying microseismic activity contained in 

record Rp^ is obtained by inverse FFT from the part Pp,„ relevant to passive 
contributions. 

rOOSn ^Updating of ratio Ep^ between the current contribution and the spectral 

model Mp.„ is performed in a relatively fast way so as to allow fast variations of the 
signal. Updating of the initial model is performed in a much slower way; in fact, the 
initial model might not be updated. 

r00521 ^The interest of these updatings is to optimize estimation of the active 

signal when its variations cannot be controlled (seasonal and daily variations), hi the 
most favourable favourable cases, they are not necessary and estimation of the active 
part' amounts to the initial model. 
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[00531 ^Another possible simplification consists in disregarding disregards the 

seismic noise at the emitted fi-equencies e mitt e d (and, if need be, at the harmonics of 
these fi-equencies). One no longer undertakes to estimat e estimation of t he active part 
and one merely cancels or interpolates the lines corresponding to the frequencies 
emitted (and, if need be, those of their harmonics) in order to obtain the passive part. 

r0054] The harmonic lines of the frequency spectrum of each signal recorded are 

first suppressed, and the amplitude and the phase are interpolated at the line values at 
the frequency of the source; the record that has not been assigned by the emissions of 
the sources possibly containing a microseismic signal is then reconstmcted by 
inversion. This implementation mode is rather suitable in the case of a limited number 
of sources and/or sources having a high signal variability with time. 



